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Quantum shape-phase transitions in finite nuclei are considered in the framework of 
the interacting boson model. Critical-point Hamiltonians for first- and second-order tran- 
sitions are identified by resolving them into intrinsic and collective parts. Suitable wave 
functions and finite-N estimates for observables at the critical-points are derived. 

> 

£3 i 1. QUANTUM PHASE TRANSITIONS 

O _ 

Quantum phase transitions are structural changes occurring at zero temperature as a 
function of a coupling constant and driven by quantum fluctuations. An example is the 
phase transition between different shapes of nuclei for which the equilibrium deformation 
plays the role of an order parameter. A proper description of such phase transitions ne- 
cessitates the inclusion of correlation effects beyond the mean-field level. In a microscopic 
approach this is done by restoring broken symmetries of the self-consistent mean-field 
solutions (HFB or HF+BCS) via angular momentum and number projection, and by con- 
figuration mixing using the generator coordinate method [ [T] • Shape-phase transitions 
in nuclei have also been studied in the geometric framework of a Bohr Hamiltonian for 
macroscopic quadrupole shapes. Analytic solutions, called E(5) [[2] and X(5) [[3], were 
obtained under certain approximations (e.g. infinite square- well potentials), and shown 
to be relevant to nuclei at the critical-points of second- and (low-barrier) first-order tran- 
sitions respectively [ H] . A key issue near criticality is to understand the modifications 
brought in by the fact that nuclei consist of a finite number of nucleons. This aspect is 
particularly important for sorting out the phase diagram and structure-evolution of stable 
and unstable nuclei, as a function of protons and neutrons. In the present contribution 
we examine this question in the framework of the interacting boson model (IBM) [5] 
which describes low-lying quadrupole collective states in nuclei in terms of a system of N 
monopole (s) and quadrupole (d) bosons representing valence nucleon pairs. This alge- 
braic model is simple yet rich enough to encompass the dynamics of quantum shape-phase 
transitions in nuclei, and illuminate the underlying finite-N structure at the critical-point. 

2. GEOMETRY AND CRITICAL-POINT HAMILTONIANS 

The starting point for analyzing shapes in the IBM is an energy surface defined by the 
expectation value of the Hamiltonian in the coherent (intrinsic) state [ El [7] 

|/3, 7 ; N) = (N\)-^(bt) N \0) , (1) 
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where tf c = (1 + /3 2 ) _1/2 [/? cos 74 + /? sin 7(4 + <#_ 2 )/y/2 + s+]. For the general IBM 
Hamiltonian with one- and two-body interactions, the energy surface takes the form 

E N (fi,j) =E + N(N - 1) [a(3 2 - 6/3 3 cos3 7 + c/3 4 ] (1 + /3 2 )" 2 . (2) 

The coefficients E Q , a, b, c involve particular linear combinations of the Hamiltonian's pa- 
rameters [[HI E] - The equilibrium shape for a given Hamiltonian is determined by the global 
minimum, (f3 e q,leq), of the energy surface, and can be spherical {J3 eq = 0) or deformed 
(Peq > 0) with prolate {b > 0), oblate {b < 0) and 7-independent (6 = 0) character. For 
(3 > the intrinsic state of Eq. ([!]) is deformed, and represents a band whose rotational 
members are obtained by projection. For example, for 7 = the band consists of states 
with good 0(3) symmetry L = 0,2,4, ... , 2N, given by 

\(3;N,L,M) ex V LM \P, 7 = 0; N) . (3) 

The projected states with fixed iV and L, involve a mixture of components \N ,n^,r, L) 
with quantum numbers related to the U(6) D U(5) D 0(5) D 0(3) chain. For [3 = 
one recovers the U(5) spherical ground state, \s ) = \N, = t = L = 0). When the 
deformed shape is 7-independent, a projection on 0(5) symmetry r is required 

\(3;N,r,L,M) oc V t ,lm\P, 1\ N) . (4) 

Stable shapes are characterized by a deep well-localized single minimum in the energy 
surface. In this case, the projection after variation procedure mentioned above provides 
an accurate approximation to members of the ground band. The dynamical symmetry 
limits of the IBM are examples of such stable structures, corresponding to a U (5) spherical 
vibrator (f3 eq = 0), SU (3) deformed-rotor [(/3 e? = a/2, 7 e(? = 0)] and an 0(6) 7-independent 
rotor (/3 eq = 1). At a phase transition the energy surface serves as a Landau's potential 
and displays a different topology. In a first-order phase transition it has two coexisting 
minima which become degenerate at the critical-point. In a second-order phase transition 
the energy surface changes continuously from one minimum to another, becoming flat- 
bottomed at the critical-point. Typical examples are shown in Fig. 1. The corresponding 
conditions on the energy surface at the critical-points are 

I s * order b 2 = 4ac, a > 0, b^O 

2 nd order a = 0, b = 0, c> . (5) 

Phase transitions for finite N can be studied in the IBM by an Hamiltonian of the 
form, Hi + gH 2 , involving terms from different dynamical symmetry chains [|7j. The 
critical value of the control parameter g c is determined by the corresponding condition in 
Eq. §5§. Several studies of this type have identified a U(5)-SU(3) first-order transition 
with an extremely low-barrier, a U(5)-0(6) second-order transition, and a 0(6)-SU(3) 
cross-over [ [TJ [10]. To consider shape-phase transitions of a general character, without 
being restricted to a particular form for the Hamiltonian, it is convenient to resolve the 
critical-Hamiltonian into intrinsic and collective parts [ [8], [9] , 

H cr i = Hi n t + H c . (6) 
The intrinsic part (H int ) is defined to have the equilibrium condensate | (3 = /3 e q,leq', N), 



Quantum shape-phase transitions in finite nuclei 



3 




Eq. ([I]), as an exact zero-energy eigenstate and to have an energy surface with the same 
shape as the critical energy surface. The collective part (H c ) is composed of kinetic terms 
which do not affect the shape of the energy surface, and can be transcribed in the form 



H c = c 3 C (3) - 6h d + c 5 Co(5) - 4n d + c 6 <%^y -5N + E 



(7) 



Here N = hd + h S) fid and h s are the total-boson, <i-boson and s-boson number operators 
respectively Cq denotes the quadratic Casimir operator of the group G as defined in [ 
[9]. The intrinsic-collective resolution constitutes an efficient method for studying shape- 
phase transitions, since the derived Hamiltonian is tailored to reproduce a given energy 
surface which, in-turn, governs the nature of the phase transition. 



3. FIRST-ORDER CRITICAL-POINT 

The energy surface at the critical-point of a first-order transition between spherical and 
prolate-deformed shapes has the following form for 7 = 

E cn ((3) = E + cN(N-l)f((3) , f((3)= /3 2 (1 + /? 2 )~ 2 (/? - /5 ) 2 . (8) 

As shown in Fig. 1(a), E cri (f3) has two degenerate minima at f3 — and (3 = /3q = 2a/b = 
b/2c > 0. The value of (3q controls the position, (3 = (3 + = (—1 + \Jl + (3q ) / (3q, and height 

h = f(P+) = (—1 + yl + Po ) 2 /4 of the barrier. In this case, the intrinsic part (H int ) has 
the equilibrium condensate | (5 — (3q, 7 = 0; N), Eq. ([T|), as an exact zero-energy eigenstate 
and has an energy surface as in Eq. ([8]) (with E = 0). H int has the form [ [TT] 

H int = h 2 P}(p )-P 2 ((3 ) , PUPo)=Pos^dl + ^7/2 (dU^)® , (9) 

with i\t(A)) = ( — l) M -P2,-/i(/3o) and h 2 = c. By construction H int is rotational- scalar and 
has the L-projected states \(3 = (3q; N, L) of Eq. ([3]) as solvable deformed eigenstates with 
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Figure 2. Left: spectrum of H int , Eq. (jSJ), with h 2 = 0.1, /3 = 1.3 and JV = 10. Right: 
the number of d bosons {n d ) probability distribution for selected eigenstates of H int . 



energy E = 0. It has also solvable spherical eigenstates: \N,rid = r = L = 0) = \s ) 
and \N,rid = r = L = 3) with energy 22 = and E = 3h 2 \J3q(N — 3) + 5] respectively. 
Accordingly, the exact spectrum of Hi nt shown in Fig. 2, displays a zero-energy deformed 
(K = 0) ground band, degenerate with a spherical (rid = 0) ground state. The remaining 
states are either predominantly spherical, or deformed states arranged in several excited 
K = bands below the 7 band. The coexistence of spherical and deformed states is 
evident in the right portion of Fig. 2, which shows the rid decomposition of wave functions 
of selected eigenstates of H int . The "deformed" states show a broad rid distribution typical 
of a deformed rotor structure. The "spherical" states show the characteristic dominance 
of single rid components that one would expect for a spherical vibrator. Table [T] shows the 
effect of different rotational terms in H c , Eq. (j7j). For a high-barrier ((3 Q = 1.3, h = 0.1), 
the calculated spectrum resembles a rigid-rotor (E ~ clnL(L + 1)) for the C3-term, a rotor 
with centrifugal stretching (E ~ a^L(L + 1) — b^[L(L + l)] 2 ) for the cs-term, and a X(5)- 
like spectrum for the C6-term. In all cases the B(E2) values are close to the rigid-rotor 
Alaga values. Insight of the underlying structure at the critical-point can be gained by 
examining the L-projected energy surface, obtained by the matrix element of the critical 



Hamiltonian © in the states ©, E { "\(3) = (/3; N, L\H cri \f3; N, L) = E^\l3) + E 



h 2 (/3-A)) 2 £ 



2,L + c 3 



+ C 6 



L(L + 1 
N(N - 1) - (1 + p 



+ c 5 



(N) 



(TV) 
2,L 



(10) 



Here d[ N ^ , S^j , and denote the expectation values in the states \/3; N, L) of ha, 
n s (h s — 1), hdifid — 1) and h s hd respectively. The mixing between the coexisting spherical 
and deformed L = states, \<f>i) = \s N ) and |0 2 ) = \/3;N,L = 0), can be studied by 
transforming to an orthonormal basis 



ri2 



rL)- 1/2 (l0 2 )-n 2 |0i)) 



(11) 
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Table 1 

Excitation energies (in units of E(2f) = 1) and B(E2) values [in units of B(E2;2t — ► 
+ = 1)] for first-order critical-point Hamiltonians [high barrier: Eqs. flS]) , (J7|) , (|2J) , low 
barrier: Eq. (|13[) ]. with N = 10. The entries in square brackets [. . .] are estimates based 
on the //-projected states, Eq. ([3]), with values of (3 determined by the global minimum 
of the respective lowest eigenvalue of the potential matrix, Eqs. ffl2l) . The rigid-rotor and 
X(5) [[3] predictions are shown for comparison. Adapted from [ [Til [T2] . 

high-barrier (/3q = 1.3, h — 0.1) low-barrier rotor X(5) 





c 3 /h 2 


= 0.05 


c 5 /h 2 


= 0.1 
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E(6f) 
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2.64 


[3.30] 




5.67 


4+ -2+ 
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1.40 
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[1.48] 


1.48 


[1.48] 
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[1.45] 
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[1.51] 


1.92 


[1.87] 
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1.37 


[1.37] 


1.37 


[1.37] 


1.44 


[1.42] 
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[1.86] 
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2.61 


0+ - 2t 


0.003 


[0.003] 


0.003 


[0.004] 


0.24 


[0.18] 


0.78 


[0.61] 




0.63 



and examining the 2x2 potential energy matrix, Kij(/3) = H cr i \^fj), which reads 
K u ((3) = E , K 12 ((3) = -c 6 (3 2 N(N-l)(l-r 2 12 )- 1 / 2 r 12 , 

K 22 (P) = [Ei N J ((3) + 2c 6 (3 2 N(N-l)r 2 12 }(l-r 2 12 r 1 + E . (12) 

The derived eigenvalues of the matrix serve as eigenpotentials, Ej~? {/3), and the corre- 
sponding eigenvectors, |$£^o)> are identified with the ground (0+) and excited (Of) L = 
states. The deformed states \j3;N,L) of Eq. with L > are identified with excited 
members of the ground-band {Lt) with energies given by E { L N \(3), Eq. (P2J). B(E2) values 
between these states can also be evaluated in closed form [ [11] [12] . The parameter (3 in 
the indicated wave functions and matrix elements is chosen at the global minimum of 
the lowest eigenvalue, E^2 ((3), of the matrix Kij(f3) (fT2l . As is evident from Tabled] 
this procedure leads to accurate finite-N estimates of observables at the critical-point. 
The characteristic spectra, discussed above, of the rotational terms in H c (0) can now be 
understood from the L-dependence of their respective contributions to E^ \{3), Eq. ( fTOl) . 
and Kij(P), Eq. (TT2|) . All rotational terms contribute diagonal L{L + l)-type splitting. 
The C6-term controls the mixing which is essential for obtaining an X(5)-like spectrum. 

The same type of analysis can be done [ [12] for a first-order critical-point with a low 
barrier. A representative Hamiltonian for this class is the critical U(5)-SU(3) Hamiltonian 

H cri = en d -KQ-Q , e = 9k{2N - 3)/4 . (13) 

Its intrinsic and collective resolution corresponds to the choice h 2 = An, P = 1/2^/2, 
c 3 = 15k/8, c 5 = — 9k/2, cq — k, Eq — —5kN in Eqs. ©,([9]). The resulting barrier is 
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extremely low, h « 10 3 . The spectrum and E2 rates, shown in Table [TJ resemble the 
X(5) predictions albeit finite-N modifications. 

4. SECOND-ORDER CRITICAL-POINT 

The energy surface for a second-order transition between spherical and 7-unstable de- 
formed shapes is independent of 7 and has the following form at the critical-point 



E cri {(3) = E + cN(N-l)g([3) , g{(3) = (3 4 (1 + (3 2 y 



(14) 



As shown in Fig 1(b), E cri (/3) exhibits a flat-bottomed behavior (~ (3 4 ) for small /3, hence 
the single minimum at (3 = is not well-localized and experiences large (3 fluctuations. 
The relevant critical Hamiltonian ([6]) has 0(5) symmetry. Its intrinsic part is given by 



Hint = t Q h d {h d - 1) 



(15) 



with to = c, and the collective part has the form given in Eq. (j7|). The states N, r, L) 
of Eq. (J4]) constitute suitable wave functions for yrast states at the critical-point, with the 
value of (3 chosen at the global minimum of the 0(5)-projected energy surface, E^2(/3) = 
{(3;N,T,L\H cri \(3;N,r,L), 



L(L + 1) - 6D[ N J + c 5 r(r + 3) - AD^ r 



+ c 6 



N{N - 1) - (1 + f3 2 ) 2 Sff 



+ E 



(16) 



Here D[^. D^f and denote the expectation values in the states \f3;N,T } L} of 

n d , fidifid — 1) and h s {h s — 1) respectively. This procedure has been tested in [ [T3] for 
a representative Hamiltonian of this class, namely, the U(5)-0(6) critical Hamiltonian, 
H cr i = ehd + APq, with e = (N — I) A. Its intrinsic-collective resolution corresponds to 
the choice t Q = A, c 3 = 0, c 5 = -A/2, c 6 = A/A, E = AN(N - l)/4. in Eqs. ©, (USD- 
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